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Abstract 

Meinardus proved a general theorem about the asymptotics of the 
number of weighted partitions, when the Dirichlet generating function 
for weights has a single pole on the positive real axis. Continuing [3], 
we derive asymptotics for the numbers of three basic types of decom- 
posable combinatorial structures (or, equivalently, ideal gas models in 
statistical mechanics) of size n, when their Dirichlet generating func- 
tions have multiple simple poles on the positive real axis. Examples 
to which our theorem applies include ones related to vector partitions 
and quantum field theory. Our asymptotic formula for the number 
of weighted partitions disproves the belief accepted in the physics lit- 
erature that the main term in the asymptotics is determined by the 
rightmost pole. 



1 Statement, history and motivation of the 
main result 

The goal of our work is the extension of Meinardus' theorem to the case when 
the Dirichlet generating function for weights has more than one simple pole 
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on the real axis. For this purpose we combine Meinardus' approach with the 
probabihstic method of Khintchine, working in the framework of the unified 
method of asymptotic enumeration of decomposable combinatorial structures 
developed in 0. 

The paper consists of two sections. In Section 1 we give a mathematical 
set up of the problem, the description of the probabilistic method used, the 
statement of our main result and a sketch of the relevant history. Section 2 
is devoted to the proofs. 

Our starting point is the following formalism. Let / be a generating 
function of a nonnegative sequence {c„, n > 0, cq = 1}: 

f{z) = Y,CnZr (1) 
n>0 

A specific feature of structures considered is that / has the following multi- 
plicative form: 

/=n^^' (2) 

k>l 

where is a generating function for some nonnegative sequence {dk{i), j > 
0, k> 1}, i.e. 

Sk{z) = J2dk{j)z'^ k>l, (3) 
j>o 

where the sequence {dk(j)} is such that the infinite product / in ([2]) has 
a nonzero radius of convergence. The above setting induces a sequence of 
probability measures n > 1 on the sequence of sets fin, n > 1 of integer 
partitions of n, such that c„ is a partition function of fin- Indeed, denoting 
by rjn = (ji, . . . , jn) G r2„ a partition of n in the sense that Ylk=i ^3k = n, 
we have from ([I]) -([3]) 

n 

Cn = ^ Jj4(jfc). 

Vn^^n k=l 

Consequently, the corresponding measures fin are defined by 

n 

f^niVn) = C'^YldkUk)- 

k=l 

The measures are called multiplicative (Vershik [11]) or Gibbsian (Pitman 
[TU]). To explain the latter name we may view /i„ as an equilibrium distri- 
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bution of a physical system, presenting /i„ in the form 

n 

KVn) = C~'^ exp ( ^ log 4(jfc)) , Vn G ^n, 
k=l 

so that log(ifc(jfc) is the total inner energy of the system at the state 

rjn. Since the above expression for the total energy does not include local 
potentials ascribed to interactions between groups of particles, the Gibbs 
measures /i„ characterize systems without interactions. 

The measures /i„ are associated with a realm of models in statistical 
mechanics and in combinatorics. The objective of the present paper is the 
asymptotics of c„, as n — oo. Our subsequent asymptotic analysis is based 
on the fundamental representation of c„ for general multiplicative measures 

Cn = e^'Ue-')F{U^ = n), n > 1, (4) 
where (5 > is a free parameter, 

n 

/n = n (5) 

k=l 

is the associated truncated generating function and, finally, 

71 

Un '■= Yk, 

k=l 

where Yk are independent integer-valued random variables distributed in ac- 
cordance with the given generating function /: 

nYk = 3k) = , J>0, k>l. (6) 

The representation (jlj) which was derived in [5] is an analogue of Khint- 
chine's representation of partition functions c„ for particular models of sta- 
tistical mechanics. In accordance with Khintchine's probabilistic method, we 
will choose the free parameter 6 = 6n to be the solution of the equation 

Ef/„, = n, n > 1, 
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which, by virtue of the above setting, can be written as 

- (\og{Ue~')))' =n, n>l, (7) 

where the derivative is taken with respect to 6. By the laws of thermody- 
namics, the function $(5) := log ^e"^/„(e~'') j is the entropy of the physical 

system considered. It is simple to show that $"(5) > 0, 5 > which tells 
us that the solution 6n of ([7]) is the point of minimum of the entropy $(5). 
This important meaning of the aforementioned choice of the free parameter 
was revealed by Khintchine ([S], Ch.6), for simple models of statistical me- 
chanics, allowing him to provide a probabilistic derivation of the second law 
of thermodynamics. 

Our study is restricted to three classic generating functions f^^\ i = 1,2,3 
associated with three types of decomposable combinatorial structures, each 
of which can be interpreted as a model of ideal gas from statistical mechanics. 
These structures are weighted partitions (Bose Einstein statistics), selections 
(Fermi-Dirac statistics) and assemblies (Maxwell-Boltzmann statistics), re- 
spectively. Defining the functions J^^^\6) = /^*^ {^~^) ; 5 > 0, for z = 1,2,3, 
we have 

^«(5) 

Clearly, by definition (|5]), the quantities fn\e~^), i = 1, 2, 3 are given by the 
expressions (15])-f ll0p . with products restricted to the range 1 < k < n. Cor- 
respondingly, the distributions for the random variables in (jl]) are of the 

following three types: Negative Binomial {bk;e~^^^ , Binomial ibk, I'^^-ik ) 



= - e-'')-'\ (8) 

k>l 



= (9) 

k>l 
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and Poisson (ft^e , and the equation ([7j) takes the forms: 



^ = ^^^^ 

fc=i 1 - e '^^^ 



Y^kb^e-''"'" = n. (13) 
fc=i 



A common feature of the three models considered is that each of them is 
given by one sequence of parameters, which is {bk > 0, k > 1}. In the 
context of statistical mechanics b^ is the weight of the energy level k, while 
in combinatorics bk is the weight prescribed to a indecomposable component 
of size k. 

Following Meinardus' approach we define two generating functions for the 
sequence {bk}'- the Dirichlet generating function D and the power generating 
function G, given by 

D{s) = J2bkk~\ s = a + it, (14) 

k=l 

oo 

G(^) = ^6fc^^ \z\ < 1. (15) 

k=l 

We assume that D and G satisfy the conditions (/) — (///) below: 

Condition (/). The Dirichlet generating function D{s), s = a + it is 
analytic in the half-plane a > p,. > and it has r > 1 simple poles at positions 
< pi < p2 < ■ ■ ■ < Pr, with positive residues Ai, A2, . . . , Ar respectively. 
Moreover, there is a constant < Cq < 1, such that the function D{s), 
s = a + it, has a meromorphic continuation to the half-plane 

n = {s:a> -Co} (16) 

on which it is analytic except for the above r simple poles. 
Condition (//). There is a constant Ci > such that 

D{s) = O {\tf') , t^oc (17) 
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uniformly for s = a + it ^ Ti. 

Condition (III). For 6 > small enough and any e > 0, 

oo 

2Y,he^''^ sm^nka) > (^1 + ^ + e) log5|, 

k=l 

v^<|«|<l/2, i = 1,2,3, 
where the constants Ai^^^ are defined by 

M^"^ = U, a 1 = 2, 

[l, if z = 3. 

Remark. The main difference between conditions (/) — (///) and the 
conditions used in the original paper of Meinardus [7] is that condition (J) al- 
lows D{s) to have multiple poles on the positive real axis. Another difference 
between our conditions and the conditions used by Meinardus is condition 
(//J), which as was shown in [3], is weaker than the corresponding condition 
in [7]. 

To formulate our main result we need some more notations. Define the 
finite set 



dkiPr -Pk): 4 e z+, 5^ 4 > 2 I n (0, pr + 1] , (18) 

fc=0 fc=o J 

where we have set po = and let Z_|. denote the set of nonnegative integers. 
Let < «! < a2 < . . . < a|x^| < Pr + 1 be all ordered numbers forming the 
set T^. Clearly, ai = 2{pr — Pr-i), if the set is not empty. We also define 
the finite set 

= t,U{p,-pfc : A; = 0,1,..., r-1}, (19) 

observing that some of the differences Pr — Pk, k = 0, ... ,r — 1 may fall 
into the set T^. We let < Ai < A2 < . . . < A|x^| be all ordered numbers 
forming the set T^. In addition to its appearance in f|T9l) . the set plays 
an auxiliary role in the proof of Proposition [1] below. 

Now that we have defined the quantities A^, we are ready to state our 
main result. 
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Theorem 1 (iMeinardus type theorem with multiple singularities 

Suppose that the conditions (I)- (III) hold. Then the following asymptotic 
formulae for Cn\ 2 = 1, 2, 3 are valid: 

1=0 1=0 s:Xs<pi 

n^oo, (20) 

where 

- { I: 

where the powers ^^-pf- do not depend on n and where the coefficients H^^\ -^l*;' ^! 
i/iai do not depend on n either, are implicitly defined in the course of the proof 
of the theorem in Section 2 below. 

A brief history and motivation. The famous Meinardus' theorem pub- 
lished in 1954, in [7], provided an asymptotic formula for Cn \ under three 
assumptions, call them (/') — (HI'), which are prototypes of our conditions 
(/) — {III). The assumption (/') differs from (/) in that it requires that D has 
only one simple positive pole, the assumption (//') is identical to (//), while 
the technical condition (///) is a weaker form of {III'). The asymptotic 
formula of [7] preceded by the great Hardy- Ramanuj an exact asymptotic ex- 
pansion (1918) for the number of integer partitions, was widely discussed in 
the literature on enumerative combinatorics and statistical mechanics. 

A formal motivation for our extension of the original Meinardus' theorem 
to the case when D has many poles comes from the basic restriction of the 
theorem, implied by its condition (/'). Due to this condition the theorem 
being valid for the weight sequences of the form 6^ = ak^~^, a, r > 0, fails 
for linear combinations of such b^s, because the induced functions D have 
many poles. 

In modern physics, the need of such extension arises in the counting of 
BPS operators, a line of research that originated in field theory in the 2000's 
(see |6]). In the language of combinatorics, problems treated in this context 
(see [H Ej) belong to weighted partitions with weights of the form bk = C^^') , 
for some / > 2. The Dirichlet generating function of such a weighted partition 
has r = I + 1 simple poles at integer points 1, 2, ...,/ + 1. 
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In the physics hterature it has been accepted the behef that in the case 
of r > 1 poles < pi,...,pr the number of weighted partitions Cn^ = 
Cn\pi, . . . , Pr) is asymptotically equivalent to the number Cn\pr) of parti- 
tions corresponding to the model with one rightmost pole pr (see (6.9) in |1] 
and (5.19) in [6]). Our asymptotic formula ( 120|) for Cn\ = 1,2,3 disproves 
this belief for all three types of the structures considered. Nevertheless, the 
basic idea behind this belief is correct, in the sense that 

log4^)(pi, . . . , Pr) ~ \og4\pr) ~ Pj^^n^, n ^ oo, (21) 

by virtue of ( 1201) . The proof of Theorem [T] shows that 

Pj*) > 0, i = 1, 2, 3. 

Finally, we note that an interesting extension of Meinardus' theorem in 
a related but a quite different direction was obtained in [8] and [9]. In both 
papers it is studied the problem of asymptotic enumeration of the number of 
unweighted partitions {bk = 1, k > 1) oi n with summands belonging to a 
given infinite set of positive integers A = {Ai, . . . , A^, . . .} and with associated 
spectral zeta function D{s) = J2i>i "^T^- the set A represents numbers 
with missing digits, so that the generating function D has simple equidistant 
poles a + 27riku, k & Z, for a fixed a; > 0, on the line 3?(s) = a for some 
a > 0, and it is analytic in the half-plane 3ft(s) > a. In [9] motivated by the 
enumeration problems in the setting of quantum mechanics, the function D 
has simple real poles at the integers n,n — 1,...,1,— 1,.... 

2 Proofs. 

The most difficult part of the proof is contained in the Subsection 2.2 in 
which we establish the asymptotic formulae for the solutions 6n\ i = 1, 2, 3 
in the case of r > 1 simple poles. 

In what follows we always suppose that conditions (/) — (///) are satisfied. 

2.1 Asymptotics of generating functions, as (5 — > 0+. 
Lemma 1 (i) As S 0^ , 

^«((5) = exp hf^S^P^ - (P)(0) log5)l(z) + M»((5; Co)j , (22) 
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where 



Po = 0, l(z) 



1, z=l; 
0, 1=2,3, 



AiT{pi)api + i), 1=1, 1=1, 

A;r(pO(l-2-^OC(A + l), ^=2, 1=1, 
AiT{pi), 1=3, 1=1, 



,r; 



D'iO), 1=1- 
L'(0)log2, i=2] 
D{0), 1=3 



and 



Co) 
Co) 
Co) 



-Co+ioD 

S-'T{s)C{s + l)D{s)ds, 

'Co—ioo 
'—Co+ioD 

rT(s)(l - 2-')C{s + l)D{s)ds, 

1 /•—Co+ioo 

— I 6~'T{s)D{s)ds, 

-Co—ioo 



1 

27ri 
1 

2Tn 



2m 



with M^''>{5;Co) = OiS^"), i = 1,2,3. 

(a) The asymptotic expressions for the derivatives 



(k) 



1,2,3, A; = 1,2,3 



are given by the formal differentiation of the logarithm of ^EEj, with 



(23) 
(24) 
(25) 



1,2,3, k = 1,2,3. 



Proof The proof of the claim {€) which we only sketch is similar to the 
proof of the first part of Lemma 2 in [3]. Following the Meinardus approach, 
we will use the fact that e~", u > 0, is the Mellin transform of the Gamma 
function: 



nv+ioo 

" = — / u-'T{s)ds, M > 0, 

^TT? Jv-ioo 



v>0. 



(26) 



Expanding log(l — e ^^) and log(l + e ^'') from (jH]) and (jO]), respectively, 
in terms of e~^'^ and then substituting ( 126|) in (E]) - (fTOl) . with v = pr + e, 
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with any e > 0, gives the desired integral representations of the functions 
log i = 1,2,3 for all 5 > 0: 

1 pe+pr+ioo 

log^«(5) = — / 6-'T{s)as + l)D{s)ds, (27) 

1 pe+pr+ioo 

logJ-(2)(5) = — / 6-^T{s){l-2-^)as + l)D{s)ds, (28) 

1 pt+pr+ioo 

log J^(=')(5) = 7^ 6^'r{s)D{s)ds. (29) 

Next, we apply the residue theorem for the above integrals, in the complex 
domain — Cq < 3ft(s) < Pr + with < Cq < 1, e > 0. By virtue of condition 
(/), the integrands in (p7|) - (!29|) have in the aforementioned domain r simple 
poles at p; > 0, / = 1, . . . , r. The corresponding residues at s = p; in the 
three cases are equal to: 

Ai5-PT{pi)aPi + l)J = l.....r, for z = 1, 

A;5-^T(p0(1-2-^0C(a + 1), / = for z = 2, 

Ai5-P^T{pi)J = l,...,r, for z = 3. 

Besides, from the Laurent expansions at s = of the Riemann Zeta 
function C(s + 1) = - + 7 + . . . and the Gamma function r(s) = ^ — 7 + . . . , 
where 7 is Euler's constant, the integrands have also a pole at s = 0, which 
is a simple one in the cases z = 2, 3 and is of a second order in the case i = 1. 
The residues at s = are equal to 

D'{0) - D{0)\ogS for i = 1, 
L'(0)log2 for i = 2, 
D{0) for i = 3. 

Finally, applying condition {II) shows that in all three cases the integrals 
of the above integrands, over the horizontal contour —Co < 3fJ(s) < e + pr, 
'^{s) = t, tend to zero, as t — j- 00, for any fixed 6 > 0. This gives the 
claimed formulae ( 122|) for z = 1,2, 3, where the remainder terms M^^^ are 
integrals of the aforementioned integrands, taken over the vertical contour 
—Co + it, —00 < t < 00. 
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The proof of (ii) proceeds as in the proof of Lemma 2 in [3]. ■ 
The asymptotics of the generating functions for weighted partitions has 
almost a century long history of its own which is enlightened in the remark- 
able monograph p. 228-229, by Korevaar, (see also [3]). 

2.2 Asymptotics of 6n- 

Preliminaries. First we note that for i = 1,2,3 equations (ITT]) - ( fT3|) have 
unique solutions 6n\ such that 6n^ — )■ as n — )■ oo. The proof of this is the 
same as the proof of part (ii) of Lemma 2 in [3]. We call any 6n, such that 

(-log/„(e-^));^^^-n^O, n^oo (30) 

an asymptotic solution of ([7]). We demonstrate that in all three cases con- 
sidered it is sufficient for (!30|) that 5„ obeys the condition 

{-\ogT{6)y^^-^^-n^ 0, n ^ cx). (31) 

In fact, by Lemma [H 

{ ~ log J^^'\S)y r^h^;^PrS~'^'\ 6^0, 1 = 1,2,3, 

from which it follows that ( 13T]) implies 

~ {h^^pr)^n"^, n ^ OO, i = 1, 2, 3. (32) 

We have, for an arbitrary multiplicative measure and all 5 > 0, 

oo 

log/„(e-^) = log J-(5) - ^ogSkiS). (33) 

k=n+l 

By using the Wiener- Ikehar a theorem as in the proof of Theorem 2 of [3], 
one may show that in our setting bk = o{k''''). (In this connection note that 
for the application of the theorem only the rightmost pole matters). Using 
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this fact, we obtain from for the case i = 1: 



oo 
k=n+l 



log si:\5, 



k=n+l 



kbke 



-kSl 



(1) 



-kS 



^(1) 

n 

oil). 



n+1 1 - 



-u5. 



(in 



(n+l)<5, 



(1) 1 - e-"' 



dw 



n — oo, 



(34) 



by fl32|) . since n6n^ — )■ oo. This together with fl33l) say that fl3Tl) imphes fl30!) 
in the case i = 1. For the cases i = 2,3 the proof of the claim is similar. 
As in [3], we will use ( 13T]) to derive asymptotic expansions for ^i*'', i = 1, 2, 3 
in our setting. It is necessary to establish that the error of approximating 
the exact solution 5„ by the asymptotic solution Sn is of order o(n~^). By the 
monotonicity in 5 > (for any fixed n > 1) of the functions log fn\e~^), i = 
1,2,3 and their derivatives, it follows from fj32|) that 



<5« n^oo, z = 1,2,3. 

Applying the Mean Value Theorem, we obtain 



(35) 



log/„(e-^" 



log/„(e-^"') 



(5»-5«)flog/, 



(36) 



where 

«„e[min(5W,5«),max(5»,5«)]. 

By virtue of the definitions ([7]), fl30l) . the left hand side of fl36|) tends to 0, as 
n — oo. It follows from the preceding arguments and Lemma [H that for all 
three cases, 



l0g/n(e^"")j ~ PriPr + l)h^\S^Y'^-' 

Combining ( 136|) with ( 137|) . gives the desired estimate 



(37) 



oin 



1,2,3. 



(38) 
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The asymptotic solution. Now our efforts will be devoted to derive the 
asymptotic expansions for ^i*-*, i = 1, 2, 3 which by will provide asymp- 
totic expansions of Sn^ up to order o{n~^). For the sake of brevity, we will 
use in the course of the proof the notations 

/i«=P;/iP, / = l,...r, ^ = 1,2,3, 

and 

_ r D(0), i=l; 
^0 - \ 0, i=2,3. 

Proposition 1 Suppose that the sequence bk > 0, k > 1 is such that the 
associated Dirichlet generating function D satisfies the conditions (/) and 
{II) of Theorem [I]. Then, the asymptotic expansions of the 5n ^ up to terms 
of order o{n~^) are given by 

1 1 

5« = ^n"^ + Kfn-l^ + o(n-^), z = 1, 2, 3, (39) 



s=l 



where K^f* and Xs do not depend on n, and the powers Xs are as defined in 



Proof By [ii) of Lemma [T] we have 

r 

- log^«(5))' = Y,hfh-'^^-' + (m«(5;Co))', 2 = 1,2,3. 



(Here and in what follows (^M'^*^(5; Co)j denotes the derivative with respect 
to 5). By (132|) . the condition (I3T1) may be rewritten as 

n(5«)^'-+i - ELo hfid^Y-"' - {^^Y^' fM»(5i*^; Co) 
A(') := 



(5«)p.+i 



o(l), n^cx), z = l,2,3, (40) 
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which is equivalent to saying that the numerator in f l40p is o{n ^). In the 
course of the proof we will set 

5« = Q» + n», z = 1,2,3, (41) 

for some qI^^ ~ (/i^^^)^^-^ ~ 6^^^ and u'n^ = o0) = o(gi*^), i = 1, 2, 3 
to be determined. We will then set Zn'^ = -%y- Since the rest of the proof 
goes the same way for all three cases considered, we suppress the index (i) 
in the forthcoming expansions. To motivate the choice of Qn^ and Un\ we 
observe that the insertion of of the form (HTl) in (HOl) . and the use the 
binomial expansion gives 



fc=0 

+ j:('"^')(^Q'n^')Zn- (42) 
m>l ^ ^ 

-ew-"(e(''';'")^")- («' 

fc=0 m>l ^ ^ 

-(m(5„;C7o)) V^J: (44) 

m>0 ^ ^ 

r 

= J]/z,gr^'=+L„, (45) 



fc=0 



where L„ denotes fl42l) + fl43|) + fl44j) . Observe that if has been chosen, 
then Ln depends on Z„ alone. 

Our plan for the remainder of the proof consists of the following three 
steps: 

Step 1. Determine a. Qn ^ Sn satisfying the condition 

r 

nQ'n'^' -J2hQ^^^~^^ = o{n~'). (46) 

k=0 

Step 2. Show that, under Qn obeying (146|) . one may choose so that 
Ln defined in satisfies L„ = o(n~^). 
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Step 1 and Step 2 construct an asymptotic solution 5n of the form (HT]) . 
By (138|) . this provides the asymptotic expansion for the solutions Sn\ i = 
1, 2, 3 up to the of order o{n~^). At this point it is convenient to introduce a 
variable 

1 

z = Zn = n '"■+1 . 

Step 3. Express the obtained asymptotic formulae for 6n\ i = 1,2,3 as 
finite linear combinations of powers z. 

Step 1. For the sake of convenience, we introduce another notation. 
Given a function Qn{z) of 2, we define 

Q{z) = z-^Qn{z). (47) 

Then P6l) can be written as 

r 

gp,.+l(^) _ ^/^^gp.-P.(^)^P.-P. _ o(n-l). (48) 
fc=0 

It follows from (gH]) that Q'''-+i(0) = /i^, while dSD implies Qp-+^{z) = 0(1). 
We will verify that fHSl) is satisfied with Q given by 

r-1 

g^'-+i(2) = Qp^+\o) + Y,hkQ''~''{0)zp^-''^ + r(2), (49) 

fc=0 

where V{z) = 0, if = 0, but where otherwise V{z) is of the form 

V^(^) = 5^5„^«'", (50) 

m=l 

where the coefficients Bm do not depend on z, and where the powers ai < 
a2 < ■ ■ ■ < comprise the set defined by (fT8l) . Denote 

r-1 
k=0 

and assume that Q{z) is defined by fj^9l) for some function ^^(2;) of the form 
fISU]) with undetermined coefficients B^- Then the condition fHS]) can be 
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written as 

r-l 



V{z) -J2hkz''^~^'{Q'''"'"'{z) - QP'-P'{0)^ = o{n-^). (51) 



k=0 

Moreover, we have 



QPr-Pk(^z)-QP'-P^{0) = 

= 0(z'"-^'-^), A; = 0,...,r- 1, (52) 

where the last asymptotics follows from the definitions of P{z), V{z). Hence 
flSTj) implies V{z) = 0{z'^^P''~P''-^)), in accordance with the fact that ai = 
2{pr — Pr-i)) provided the set is not empty. 

It follows from the definition of the set that if we can find V{z) of the 
form fISUl) such that 

Bm - [^"'"] ( Yl hkz"^"'' {Q"^-'' {z) - QP^-P^ (0)) ) = 0' 

\ k=0 ) 

m = 1, . . . , |T,.|, (53) 

where the second term in flS51) is the coefficient of z"™ in the expression in 
parentheses, then we will have found the required Vi^z^ satisfying ( 15T|) . Using 
( 152|) . the equation (153|1 with m = 1, ai = 2(pr — Pr-i) can be written as 

2 (Pr-Pr.-l)g^'"'''^-^HO) 



Bx = K_, — . 

Pr + 1 

The crucial fact is that, by virtue of (15^ . the coefficient of 2;°™ in (1551) does 
not depend on Bg for s > m, so that 

Bm + qm{Bi, . . . , Bm~i) = 0, m = 2, . . . , \Tr\, 

where are polynomials of Bi, . . . , Bm~i- Thus, the Bm can be defined 
recursively from (1531) . resulting in a function V{z) of the form (l5Up satisfying 

Finally, with this choice of V{z), we obtain from fl52|) 

/ P(z) + V(z)\^ 

g(^) = Q(0) 1 + ^ ' ^ 
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Recalling ( H7|) . this gives 

Q^^K^M^mim)^^, (54) 

which, by our construction, satisfies (l46l) . 

Step 2. We first need several definitions. Under Qn given by ( |54l) . define 

C.(m) = f + ^\ uQ^^^^' ~ f ^^\k^Q, m > 1, (55) 
\ m J \ m J 

D„(m) = 5^/iJ^'^^^Mgf;-^'»=0(g^'-^'-i), m>l, n^oo, (56) 

fc=0 ^ ^ 

and 

En{m) = (^^'^ (m(5„;Co))'q^'-+\ m > 0, n ^ cx). (57) 



In accordance with fH5l) . 

i^n = $^(Cn(m) - Z}„(m) - K(m))Z™ - K(0). (58) 



m>l 



We will analyze the factor yM{Sn', C'o) j in fl57|) . Suppose that 

Z„ = O (Q^«+'"-) . (59) 

Evidently, 

1 p — Co+ico 

(M(5;Co))' = — / r^-i/i(.)d., 

where h{s) is determined by the appropriate choice of ( l23l) - ( 1251) . Under 
= Qn(l + ^n), we havc 

(m(5„;Co))' = Q^«"i(l + Z„)^"-'( rgr"(l + ^n)"'"M-Co + ^r))zrfr 
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For given < Cq < 1, Pr > as in the statement of our Theorem 1, 
define 

mo = min{m > 2 : (Co + Pr)m > pr + 1}, (60) 

so that = o(n~i). 
Next we have 

(I \ mo-2 
^ ' 1=0 

m > 0, n — )■ oo, 

where 

Vi{n) = j Q;;^(^'^^^h{-Co + iT)idT = 0{l), />0, n->cx) (61) 
are known constants. Now the expression for En{m) can be written as 

mo— 2 
1=0 

where we denoted 

7z(n,m)= (^^^ + ^^ J2 (^'p ^^Vk{n). 

Consequently, 

mo— 2 

^E„(m)Z™ = J]E„(m)Zr + o(n-i) 

m>0 m=0 

mo— 2 

= Qr^"" E E 7z(r^,m)) +o(n-^) 

fc=0 m+l=k 
mo— 2 

= Q^^^""" + (62) 

fc=0 

where 

lk{n) = ^ Tzl'^,"^)- 
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Set 

Wn{m) = Cnim)-D^{m), m > 1, (63) 
Substituting (1^ and into (IHSj) produces 

mo — 1 mo— 2 

i^n = £ W^n(m)Z- - ^ 7,(n)Z„^ + o(n-i). (64) 

m=l fc=0 

Our objective will be to show that the condition 

Ln = o(n-i) (65) 

is satisfied by Z„ which, in agreement with the assumption (1590 . is a poly- 
nomial in Q'^^^^ of degree (mo — 1) : 

mo-l 

Zn=Y. Pk{n)Q':i^^^+^^\ (66) 
fc=i 

with coefficients /3fc(n) = 0(1) that are determined recursively from (!65l) . In 
fact, under the above form of Z„, the main term of L„ becomes a polynomial 
in Q^^'"° , so that the condition is satisfied if 

[Q^(^^+^«)]L„ = 0, A; = 1,..., mo-l, (67) 

where [g^^''''+^''^]L„ denotes the coefficient of Qn^""-^^''^ in From ([671), 
(1M|) we derive 

[g(f+^«)]L„ = W^„(l)/3i(n) - 7o(n) = 0, 
[g2(p.+Co)]^^ = Wn{l)P2{n) + W„(2)/3i2(n) - 7i(n)/3i(r2) = 0, . . . , 

mo-l 

[QKf'r+O^L^ = [g^(''^+^«)]( ^ l^„(m)Zr))- 

m=l 

mo— 2 

^g(fc_i)(,,+C7o)]( ^ 7,(n)Z9 =/3fc(n)iy„(l) + i?,_i(/3i(n),.../3fc_i(n)), 

A;=0 

/c = 1, . . . mo — 1, 

where Rk-i is a polynomial in /3i(r;,), . . . Pk_i{n). By induction on fc, we have 
/3/c('^) = 0(1) for all /c G [l,mo — 1]. Since Rk-i does not depend on /3fc(n) 
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and since Wn{l) ~ (pr + l)^r 7^ as n — )■ oo, our claim is justified. Thus, 
with Zn given by 



I 7mo 



m>mo " 



by virtue of fl60|l . (|63|1 and because Z„ — )■ 0. Consequently, we have m„ = 
^nQn = 0(Q^o+'''-+^) = o(n-^), as was desired. 

Step 3 The first two steps, combined with (l38l) produce the following 
asymptotic formulae for the solutions 5^*^ : 



Pr + l / _l^ 

+ o(n , 



^ hr ^ 

(68) 

where Vi^z) is given by (1501) . We now observe the following fact stemming 
from the definition of the set T^. If A^, G T,. are such that A^ + A^ < Pr + 1, 
then As + Ai e as well. Hence, the binomial expansion of the right hand 
side of fl68|) gives the claimed formulae fl39l) .B 

The following particular cases of Proposition 1 are of importance. 

Corollaries We recall that the set is not empty if and only if cti = 

2(Pr - Pr-l) < Pr + 1- 

[%) The case when p^ > 2pr-i: Though in this case the set T,. is not 
necessary empty, (l68l) conforms to 

5« = ih^V^n^ + ^^^-"^^^^^-^^"T"""^ + o{n\ (69) 

Pr ~l~ 1 

For the proof recall that if t,. ^ 0, then V{z) ~ Biz^^P'-p^-''^), so that in 
the expansion fl68l) . n'pr+i ^(z) = o{n~^). 

(ii) The case of a single simple pole: r = 1, pi > 0. Since po = 0, 
the condition of Corollary (i) holds and, therefore 

t{i) -1 

5» = (/i«)^n-^ + + o(n-^), (70) 

Pi + 1 

which recovers the asymptotic formulae for 6n^ in [3] . In this regard we note 
that in [3] a detailed form of the term o(n~^) was derived. The latter is not 
needed for our subsequent study. 
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{iii) The case of equidistant simple poles: pi = la, I = 0, . . . ,r, a > 
is a given number. In this case the condition of Corollary (i) holds only if 
r = 1. For r > 2, the set consists of multiples of a > 0: 

Tr = {p = aM : M eZ+,M <r + [a~^J } . 

The simple structure of allows to write the expression (1681) for in the 
following form 

1=0 

where V'o*^ = 1; "^ = 15 2, 3. We will demonstrate that the coefficients z/^j'*^ are 
determined recursively from ( H6|) . In view of the above representation of 
the condition ( 146|) can be written as 

which is equivalent to 

e=o k=0 £=0 

s = 0, ...,r + [a-^\. 

This gives 

(ar + + qs{4'\^|J^\ . . . , V'Si) = 0, s = 1, . . . , r + 1, (71) 

where Qs is a polynomial of ?/'|'\ . . . , ipg^li only. ■ 

(it;) The case where V{z) influences the main term. In all of the 
previous corollaries the remainder V{z) does not influence the main terms, 
i.e. terms of order > n~^, in the asymptotics of . We give now an example, 
where the impact of V{z) is seen. Suppose r = 2 and 2(p2 — Pi) < P2- Then 
tti = 2(p2 — Pi) G T^, and moreover, zV{z) ~ zBiz°'^ ^ n^^, as n — )• oo. 



21 



2.3 Completion of the proof 

It is left to use Lemma [1] and Proposition [1] to find the asymptotics of the 
logarithms of each of the three factors in the representation (jl]) when 6 = 6n\ 
i = 1,2,3. 

n5» = + Kf^^f^ + e«, ^ = 1, 2, 3, (72) 

s:Xs<pr 

where G T^, e^n 0. 

{ii) Firstly, by the argument similar to the one for the proof of ( l34l) we 
conclude that 

log/»(e-^"') = log^«(5«) + e«, z = 1,2,3. (73) 
Next, for / = 0,l,...,r 

s:Xs<pi 

where en\l) = o (1), z = 1, 2, 3, / = 1, 2, . . . r and where the coefficients K^^j 
are obtained from the binomial expansion for (^n"*) ^\ based on ( 139|) and the 
definition f fTOj) of the set T^. Consequently, substituting 6 = 6n^ into f l22|) 
gives 

iog/i^He-^^'^)=x:/^r^(/^«)^n^+E/^F E 

(^\ogn-^\ogh^;Ali^) + e^\ . = 1,2,3. (74) 

(in) It turns out that Theorem 4 of [3], which is a local limit theorem, 
may be adapted to the situation here. We will explain in the proof below 
that the reason for it is the fact (l2Tl) . which we pointed in our historical 
remark. 

Theorem 2 (Xocal limit theorem^. 

Let the random variable Un be defined as in Section 1. Then 
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P {U^^ = n) 



27rVar(ui'^) \/2ttK. 



1_ ^^^)y+P,./2 



1 , ^ , . 2+Pr 2 + pr 

^ (/^W)2(p.+i)y^ 2U^^ n ^ cx), z = l,2,3, (75) 



wi/i constants defined by 



Ki'^ =Arr{pr + 2)aPr + l), 



Kf = Ar{l - 2-P^)r{pr + 2)C(p. + 1) 



and 



/s:f) = A,r(p, + 2). 



Proof. We will sketch the proof that follows the pattern in [3]. Denoting 



Ua) = Eie 



^2-KiciUr, 



a G 



the characteristic function of the random variable Un and setting 

«o = (5i^^)'^logn, 



we write 



where 



and 



1/2 

¥{Un = n)= I 0„(a)e-2™"cia = Ji + /2, 
-1/2 



-QO 



-ao /•1/2 



1/2 



Defining i?„ and T„ by 

52 = ( log /„ (e"^) ) and T„ = - ( log /„ (e"^) 



5=5: 



(i) 



(76) 
(77) 

(78) 
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for n fixed we have the expansion 



0„(a)e-2™" = exp (27ria(Et/„ - n) - 2Ti'^a^Bl + 0(a^)r„) 

= eyi^ {-2Ti'^a^Bl + 0{a^)Tn), a ^ 0. (79) 

By virtue of f l22p and f l73p we derive from fl78p that the main terms in the 
asymptotics for B"^ and T„ depend on the rightmost pole only: 



■2 ) 

three cases, 



where K^'',K^'' are as in the statement of the theorem. Therefore, in all 



Consequently, 



B'^al — oo, T„aQ — )■ 0, n — )■ oo. 



Ii ~ — , n — )■ oo 



and it is left to show that I2 = o(/i), n — )■ 00. Taking into account that for 
all three models ao = o{^/6^), we then follow [2] splitting the range [oq, 1/2] 
of the integral I2 into three subintervals [oq, Sn], [Sn, v^]? and 1/2]. 
The proof of Lemma 3 in [3] shows that 

< (l+£n)exp ( -^5:^^ ) , aeR, 



where — 0, A^^*-*, z = 1,2,3, are positive constants defined in condition 
{in') of Theorem 3 in [3], and 

00 

Vn{a) := ^e"^'''"sin2(7raA;). 

k=l 

In the third subinterval, the condition (///) holds, which allows to derive, 
in the same way as in [3], that the corresponding part of the integral I2 is 

O ((5^^^) = o(/i). 

Regarding the first subinterval, we are able to derive the desired bound, 
following the scheme in |2], modified so as to match our setting for the se- 
quence {bk}- The subsequent analysis uses an inequality (3.70) from [2], 



24 



which we recall in f l8Up below. Let [x] and {x} denote, respectively, the inte- 
ger and fractional parts of a real number x, and let denote the distance 
from X to the nearest integer, so 

{x} if {x} < 1/2; 
1 - {x} if {x} > 1/2. 

We then have 

sin^(7rx) > 4 || x ||^ . (80) 
For a G [ao; ^n], we get the estimate, as n — > oo: 

n 

Vn{a) = J2bke-''"\\ajf 

k=l 

l/(2(5„) 

k=l 
l/(2<5„) 

k=l 

^ «oCi-^(25n)"^^"^ (81) 
= C2 log^ n. 



> ale Yl ^'^'i 



2P-+2(p^ + 2)' 

In order to show the asymptotic equivalence flSTl) . we used the Wiener- Ikehar a 
theorem (see p. 122), which for {6^} in our setting reads as follows: 

n 

fefc/c"^"'''' ~ Aril, n — )> 00. 

fe=l 

By writing ^22=1 ^kk^ = X]fc=i(^fc^^~'''^)(^'''^^*~^)' ^"^^ s > 0, and using summation- 
by-parts, we obtain 

n . 

\^ b^k' —nP^+', n^oo, (82) 

which, with s = 2, implies fIST]) . 
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For a G [5„,, define, as in [2], 



Q{a) 



1 < A; < n : J + 1/4 < aA; < J + 3/4, j = 0, 1, 



Aan — 3 



where Qj(Q!) denotes the set of all integers k G [^|^, ^4^]- Observe that for 
any a G [5n, and j > 0, the set Qji^ot) is not empty, since in this case 



4i+3 4i+l 



4q 



4a 



> 1. Now the aforementioned inequality from [2] yields 



VJa) > 



~k5„ 



(83) 



Next, using (152]) with s = and the fact that 5„a ^ < 1, a G we 
estimate the inner sum in (l83l) . for j > : 



5^ &fce" 



> 



r 4j+3 i 
Pr 





r4j + 31 


Pr 


r4j + ii 




4a 




4a 



Pr 



+ max{l,j'''}o((l/a)'''-) 



> ::!:(4a)-'''-e-(^+=^/^) ((4j + 3)^'- - (4j + l)"-- + o(j^'-)) 

Pr 

> C5-^^/2e-^' ((4j + 3)''^ - (4j + 1)''^ + o{r)) , a G v^], 

where we denoted C = ^e~^/'^A~'^'' . Now continuing (l83ll . we get with the 
help of the Euler integral test 

r Aa.n — Z i 



-1 L 4 J 

> Y^C^^n""^/' e-^((4j + 3r-(4j + l)^^+o(/0) 

i=o 

1 /""^ 

= Ci^;:^'-/^ Ci > 0, a G [5n, n^oo. 

From the preceding analysis it is easily seen that e-^"(") = o(/„), 
for all a G [ao, v^]- ■ 
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Finally, to completely account for the influence of all r poles p, . . . , pi, we 
present the sum of the expressions (172]) , (17111 obtained for the flrst two factors 
(i) and [ii), in the following form: 

n5«+log/„(e-^^^') = ±Pl'^n^^+±hP ^1 

1=0 1=0 s:Xs<Pi 

+(m,„g„_m,og,-4..)i(,)+4.). 

V Pr + 1 Pr+l ^ 

where denotes the resulting coefficient of np^. 
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